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Abstract 

Studied  and  developed  is  a new,  relatively  simple,  rational,  and  viable 
theory  of  phenomenological  thermal  turbulence  occurring  in  incident  unstably 
stratified  fluid  flows.  The  subject  of  stratified  flow  has  consisted  largely 
of  the  study  of  stably  stratified  incident  flows  and  of  unstable  thermal  con- 
vection of  a fluid  heated  from  below  with  no  incident  flow.  The  new  turbulence 
theory  is  intermediate  in  complexity  between  the  venerable  but  generally  inappli- 
cable turbulent  eddy-transport  theory  and  the  very  complex  generic  theory  known 
as  "turbulence  modeling"  enjoying  much  current  attention.  We  call  our  method 
"a  limit  theory  of  phenomenological  thermal  turbulence  modeling";  based  on  a 
system  of  plausibly  large  and  small  parameters,  it  is  more  suitable  in  many 
applications  where  the  said  parameters  behave  as  required;  moreover,  the  theory 
lends  itself  to  systematic  improvement.  The  work  proposed  is  preceded  by  experi- 
mental probing  in  a low-speed  wind  tunnel  and  use  is  always  made  of  the  general 
equations  as  a continuing  guide  to  the  specializations  made  in  our  specific  for- 
mulation that  sidesteps  the  attendant  major  computation  in  the  general  unstable, 
turbulence  formulation.  The  research  consists  of  a study  of  turbulent  cellular  flow 
due  to  an  unstable  temperature  stratification  underneath  a flow  with  a stable 
temperature  stratification.  These  two  layers  are  separated  by  an  interface  whose 
detailed  structure  is  a pausal  boundary  layer.  The  turbulent  thermal  convection 
below  buffets  the  upper  layer  by  means  of  "penetrative  convection"  and  heats  the 
upper  layer  at  the  interface  because  the  turbulent  heat  flux  mechanism  must  switch 
off.  The  penetrative  convection  and  heat  switch  must  amount  to  a "heat  engine" 
whereby  the  rate  of  work  of  thermal  origin  is  done  on  the  bottom  of  the  upper  layer. 

This  power  source  can  then  produce  internal,  surface,  and  Kelvin-Helmholtz  waves. 
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1.  Introduction 

Studied  is  a i-platively  simple,  rational,  and  viable  theory  of  thermally 
turbulent  incident  stratified  flow,  dominantly  controlled  by  the  Rayleigh 
number  R (although  other  parameters  are  important),  in  contrast  to  the 

gI 

classical  inertially  turbulent  flow,  controlled  by  the  Reynolds  number  R^. 

The  application  given  of  the  theory  is  kept  reasonably  specific  in  order  to 
avoid  the  pretentiousness  of  studying  "the  turbulence  problem."  The  work  is 
accompanied  by  prior  experimental  probing  in  a low-speed  wind  tunnel  (see  Appendix) 
in  order  to  provide  satisfaction  that  the  theoretical  simplifications  presented 
here  are  commonly  grounded  in  laboratory  fact.  Use  is  always  made  of  the 
general  equations,  within  the  Boussinesq  approximation,  as  a continuing  guide  to 
the  specializations  made  in  the  specific  formulation  that  sidesteps  the  attendant 
major  computation  identified  with  turbulence  modeling  as  generally  practiced. 

Our  prior  theoretical  and  experimental  probing  suggests  notable  promise  for 
a new  phenomenological  turbulence  theory  intermediate  in  complexity  between  the 
venerable  but  generally  inapplicable  turbulent  eddy-transport  theory  ("K-theory") 
and  the  generic  theory  known  as  "turbulence  modeling"  enjoying  much  current 
attention.  We  call  our  method  "a  limit  theory  of  phenomenological  thermal 
turbulence  modeling";  based  on  a system  of  plausibly  large  and  small  parameters, 
it  is  more  suitable  in  many  applications  where  the  said  parameters  behave  as 
required;  moreover,  the  theory  is  subject  to  systematic  (i.e.  rational)  improve- 
ment. The  turbulence  modeling  schemes  are  very  complicated  but  are  still  flawed 
by  the  weak  link  of  any  of  the  quite  crude  and  arbitrary  closure  relations. 

Thus  we  develop  a new  theory  of  phenomenological  thermal  turbulence  occurring 
in  incident  unstably  stratified  fluid  flows;  it  contains  both  new  ideas  and  new 
methods.  If  the  unstable  flow  lies  under  a stable  flow,  a heat-switch  mechanism 
can  produce  a heat  engine  near  the  interface;  we  show  this  result  explicitly. 
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Previously  the  subject  of  stratified  flow  has  consisted  largely  of  the  study  of 
stably  stratified  incident  flows,  with  consequent  wave  phenomena,  and  of  unstable 
thermal  convection  with  no  oncoming  stream.  The  study  of  unstably  stratified 
incident  flow  has  not  received  commensurate  attention  (although  its  importance 
is  adequately  known)  presumably  because  one  must  cope  simultaneously  (as  we  do 
here)  with  the  cells,  turbulence,  and  the  oncoming  stream.  Thus  the  objective 
of  the  proposed  work  is  to  do  just  that.  We  provide  these  new  ideas  and  new 
methods  to  deal  with  a little  explored  area  to  develop  significant  new  results. 

There  are  many  areas  of  relevance  in  geophysical  fluid  dynamics.  For 
example,  the  unstable  turbulent  "mixed  layer"  on  the  ocean  surface  largely 
drives  the  layers  of  abyssal  sea  water  beneath.  Further,  in  addition  to  the 
tides,  the  water  above  the  continental  shelf  is  well  stirred  by  the  surface 
wind  stress  and  driven  by  the  esturine  currents.  Moreover,  these  results  have 
application  to  improved  technology  for  off-shore  oil  drilling  on  the  continental 
shelf;  near  the  surface,  near  the  bottom,  and  in  between.  Submersible  objects 
are  anticipated  at  all  three  levels  — staged  or  unstaged.  Thus  the  results  can 
provide  aid  with  respect  to  fossil  fuel  recovery.  In  the  atmosphere,  the 
surface  wind  layer  of  the  planetary  boundary  layer  is  often  unstable  and 
turbulent  and  this  fact  is  relevant  to  the  energy  and  pollution  problems  plaguing 
modern  society.  Further,  in  the  atmosphere  there  are  the  flow  complexities,  in 
the  vertical,  of  frontal  systems,  and  also  the  convective  flow  field  associated 
with  clouds  outside  the  visible  cloud  structure:  these  two  problems  are  relevant 
to  weather  forecasting. 

Conerning  the  theory,  it  is  appropriate  to  begin  with  the  original  work 
(Goldstein  1938  and  Schichting  1968)  in  turbulence  wherein  the  density  p is 
constant  and  the  temperature  T is  absent.  For  the  case  of  incident  flow,  one 
decomposes  the  pressure  and  velocity  variables  into  a time  mean  plus  an  instanta- 
neous fluctuation  thus: 
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p = p + p*  and  u.  = u.  + u!,  where  (i  = 1,2,3) 

f r r 1X1 

in  cartesian  tensor  analysis.  One  then  forms  the  time  mean  of  the  continuity 
and  momentum  equations,  the  inertia  terms  of  which  contribute  the  extra  six 
turbulent  Reynolds*  stresses 

s -p  uTuT  ; (i,j  = 1,2,3) 

These  stresses  can  be  evaluated  in  terms  of  the  local  mean  flow  by  eddy  trans- 
port considerations,  viz.,  the  Prandtl  mixing  length  theory,  von  Karman's 
hypothesis,  or  the  vorticity  transfer  theories  of  Taylor  and  Goldstein.  However, 
to  provide  instead  for  the  fact  that  the  turbulence  generally  has  dynamics  of 
its  own,  one  multiplies  the  momentum  equations  by  Uj^  and  then  forms  the  time 
mean  to  obtain  the  six  fluctuation-moment  equations;  these  equations  now  govern 
the  Reynolds*  stresses  but  contain  other  complex  mean  correlations  in  addition; 
the  problem  of  relating  these  extra  correlations  to  the  existing  mean  variables 
above,  in  order  to  close  the  system  mathematically,  is  the  original  and  perplexing 
closure  problem  by  "turbulence  modeling,"  say,  for  which  no  fully  satisfying 
arguments  exist.  These  incident  flow  considerations  of  turbulence  are  entirely 
inertial,  whereas  the  present  research  emphasizes  thermal  turbulent  correlations. 

If  we  pass  to  a flow  initially  at  rest  (no  incident  flow)  but  heated  from 
below  at  sufficient  R , then  gravitational  instability  overcomes  molecular 
transport  and  cellular  motion  ensues.  In  its  linear  version  we  have  the  famous 
B4nard  problem  (Chandrasekar  1961);  as  R is  continuously  increased,  the  flow 
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becomes  increasingly  nonlinear  (and  very  sensitive  to  Prandtl  number  P^  unless 
P^  is  very  small)  and  finally  turbulent  (Turner  1973).  The  proposed  problem  has 
antecedents  in  this  problem  of  thermal  convection. 

Current  approaches  to  the  problem  of  turbulence  in  an  incident  flow  repre- 
senting research  in  progress  by  turbulence  modeling  are  prominently  represented 
by  the  works  of  Lumley  and  Khajeh-Nouri  (1973),  Launder  and  Spalding  (1972), 


Donaldson  (1972a),  and  Daly  and  Harlow  (1970) j Deardorff  (1973)  employs  direct 
numerical  simulation  of  turbulence  but  with  turbulence  modeling  on  the  sub-grid 
scale,  where  the  grid  scale  is  chosen  no  smaller  than  practicable.  Other  cxirrent 
research  in  progress  concerning  incident  turbulence  flow  includes  direct  numeri- 
cal solution  of  the  Navier-Stokes  equations  by  Orszag  (1972). 

Current  research  in  progress  in  stably  or  unstably  stratified  incident 
turbulent  flow  conspicuously  includes  the  studies  of  Lumley  (1973,  1975), 
Donaldson  (1972b,  1973),  Tennekes  (1973a, 1973b) , and  Wyngaard  (1973).  These  works 
make  use  of  turbulence  modeling  and/or  similarity  laws.  Background  reference 
books  are  those  by  0.  M.  Phillips  (1969),  P.  Morel  (ed.)  (1973),  and  especially 
D.  A.  Haugen  (ed.)  (1973). 

The  turbulence  modeling  in  the  fluctuation -moment  equations  of  the  extra, 
complex,  mean  turbulent  correlations  is  performed  in  order  to  close  mathemati- 
cally the  system  of  relations  in  use.  All  of  the  known  turbulence  modelings  are 
rather  arbitrary  and  questionable,  and  certainly  very  crude  when  compared  to  the 
closure  relations  of  kinetic  theory  or  even  radiative  transfer  by  differential 
approximation.  Whereas  there  is  merit  in  the  argument  (Lumley  and  Khajeh-Nouri 
1973)  that  the  crude  third-order  closure  schemes  are  attenuated  in  their 
influence  on  the  first-order  equations  of  mean  motion,  the  closure  schemes 
still  serve  to  call  in  question  the  common-sense  worth  of  the  extremely  extensive 
calculations  with  20  to  30  simultaneous  partial  differential  equations  that 
follow;  moreover,  the  computer  and  fluid-dynamic ist  team  can  raise  a serious 
question  concerning  which  one  is  slave  and  which  one  master? 

2 , The  General  Problem 

In  Sects.  2.1  and  2.2  we  present  the  general  problem  as  it  might  be  formu- 
lated by  workers  with  a hyperthetically  large  preceding  literature  and  very 
extensive  computing  and  experimental  capabilities.  In  Sects,  2.3  and  2.4  we 
present  what  we  shall  call  the  "heat  switch  mechanism"  and  the  "heat  engine," 
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which  we  show  occur  when  the  gravitationally  unstable  fluid  is  bounded  by  stable 
fluid  above.  We  later  obtain  the  heat  engine  operation  explicitly. 

2.1  Physical  Statement  of  General  Problem 

We  consider  a moderately  unstable  stratified  shear  flow  furnished  with  a 
lid  and  floor.  In  the  atmosphere  or  ocean  (salinity  aside),  we  expect  that 
this  flow  will  not  only  be  cellular  but  also  turbulent  for  thermal  buoyancy 
reasons  dominantly  controlled  by  R . 

From  the  geophysical  point  of  view,  a relevant  unstably  stratified  shear 
flow  is  the  one  shown  in  Fig.  1 where  the  interacting  lid  is  a stably  stratified 
shear  flow  above  and  the  floor  is  a solid  boundary  endowed  with  temperature 
or  heat  transfer  rate  q.  that  drive  or  do  not  drive  the  initial  flow;  the  non- 
monotonic  temperature  profile  TCx^)  and  velocity  profile  U(Xg),  together  with 
a turbulence  field  are  introduced  initially  at  x^  = 0.  This  problem  is  clearly 
difficult  but  the  unstable  layer  can  represent  the  following:  (i)  in  the  atmosphere, 
an  unstable  earth  surface  layer  interacting  with  a stable  upper  planetary  boundary 
layer  and  possible  producing  internal  waves,  Kelvin-Helmholtz  waves  or  just  a 
general  lift  of  the  interface,  and  (ii)  the  upper  ocean  (inverted),  where  the 
mixed  layer  (with  incompletely  modeled  surface  wind  effects)  interacts  with  the 
stably  stratified  fluid  below  and  can  produce  internal  waves  and  Kelvin-Helmholtz 
waves  again.  This  problem  introduces  the  notable  physical  phenomena  discussed  in 
Sects.  2.3  and  2.4  below. 


2.2  Formulation  of  the  General  Problem 


Let  be  the  velocity,  p be  the  density,  p be  the  pressure,  T be 

the  temperature,  B be  the  thermal  coefficient  of  expansion,  u be  the  viscosity, 
and  ( )^  denote  the  constant  common  reference  state.  The  strain  rate  is  given 
by  e^j  H + Uj,^.  In  accord  with  a tenent  of  the  Boussinesq  approximation, 

we  split  out  the  dominant  part  of  the  body  force  in  the  momentum  equation  M by 
introducing  the  neutral  reference  hydrostatic  equilibrium 
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- P^g63i  = P,i  ; (i  = 1,2,3), 

where  the  dominant  P 'v<  is  a reference  hydrostatic  pressure. 
We  introduce  the  departure  pressure  p by 

P = P + P 

and  the  departure  temperature  T by 


T = T + T, 
o ’ 


(2.2.1) 


(2.2.2) 


(2.2.3) 


where  the  tilde  C)  is  omitted  in  turbulence  equations.  The  thermal  equation 
of  state  in  the  Boussinesq  approximation  is 


p = p (1  - B T). 
o o 


(2.2.4) 


The  general  departure  equations  of  continuity  C,  momentum  M,  and 
energy  £ are  then,  for  constant  properties, 


and 


C:  u.,.  = 0,  (2.2.5) 

J J 

8:  + ^PoVj^’j  = Po^STdg.  - P,.  + p^u.,..  ; (i,j  = 1,2,3), 

(2.2.6) 

(p  c $),.  + (p  c Tu.),,  = - F.,.  + k f,..  + e..e..,  (2.2.7) 

= '^o  po  ’t  "^o  po  ] 3*3  o *33  2 *^o  13  13’ 


where  c^  is  the  specific  heat,  k is  the  thermal  conductivity,  and  is 

the  radiative  flux  vector.  These  equations  are  identified  with  unsteady  thermal 
convection. 

By  writing  u^  = u^  + u^,  T = T + T*,  etc.  and  taking  the  time  mean  of 
the  above  equations,  we  obtain  the  general  mean  turbulent  departure  equations. 
These  are 


C: 


u. , . = 0, 
3 3 


(2.2.8) 


and 


Po^®^^3i  - P’i  ^ Vi’jj  ■ 

; (i,3  = 1,2,3),  (2.2.9) 


(2.2.10) 


E:  (p  c T).  + (p  c Tu.),.  = - F.,.  + k T...  -(pc  T'u!),. 

= o po  ’t  o po  3 *3  3*3  o *33  *^0  po  3 ’3 


+ ^ye,.e..+—pe!.e!.  , 
2 o 13  13  2 o 13  13  ’ 


where  the  "extra"  terms  are  the  turbulent  correlations 


Z . . = - p u'.u! 

013 


Q.  E p c T'ul 
3 o po  3 


^ y e! .e! . = 


2 o ij  ij 


turbulent  stress  tensor  for  constant 
density  flow, 

turbulent  heat-flux  vector  for  constant 
density  flow,  and 

viscous  turbulent  dissipation  per  unit 
volume  for  constant  density  flow. 


It  is  observed  that  ten  temporal  correlations  stand  unevaluated  and  these  are 
not  taken  herein  to  depend  on  or  track  the  local  mean  flow  field;  the  problem 
is  then  very  formidable  and  not  herewith  closed  mathematically. 

We  can  form  the  general  mean  fluctuation  departure  mechanical  energy 
relation  by  multiplying  (2.2.6)  by  Uj^;  by  writing  the  variables  as  u^  = u^  + u 
T = f + T*,  etc.;  and  by  taking  the  time  mean  of  the  resulting  equation.  Then 
integrate  over  the  finite  volume  t bounded  by  the  fixed  surface  A to  obtain 
finally 


o (-  p'u!  + y e!.u!)n.dA  - p u!u!u!n.dA  = 
Ja  3 o 13  1 3 2 J/o  1133 


f »o  it  - 1 * I 


(2.2.11) 


e ! .e ! .dt  - — 1 

, 13  13  2 


Z . .e. .dt. 

,13  13 


If  is  the  departure  traction  on  A,  then  the  first  integral  on  the  left  can 

be  shown  to  be  9 T!u!dA,  where  T!  = a!.n..  In  stratified  flow  a key  point  is 
Ta  1 1 1 13  3 ^ ^ 


that  the  departure  body  force  per  unit  volume  (5  buoyancy  force)  in  (2.2.6)  is 


nonconservative;  key  points  in  the  heat-switch  mechanism  presented  next  are  that 
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T'u'  >0  in  the  unstable  flow  and  that  the  second  integral  on  the  right  makes 
o 

a large  contribution  in  the  unstable  flow. 

2.3  Heat-Switch  Mechanism 

Consider,  at  time  t = 0,  the  initial  vertical  mean  turbulent  profiles  of 
temperature  T(Xg)  and  horizontal  stream  shown  in  Fig.  2 (derived  from 

Fig.  1).  A pause  layer  is  shown  and  the  temperature  increases  from  this  value 
in  both  the  upper  layer,  denoted  by  subscript  u,  and  the  lower  layer,  denoted 
by  subscript  1.  Then  the  upper  layer  of  fluid  is  gravitationally  stable  because 
3f  /3x  > 0 and  the  lower  layer  is  gravitationally  unstable  because  3T  /3x  < 0; 

these  two  segments  of  the  initial  temperature  profile  have  no  special  shape, 
e.g.  they  are  not  straight  lines.  The  initial  horizontal  stream  profile  has  a 
shape  that  is  not  of  special  importance  to  the  present  consideration  excepting 
that  the  profile  in  the  lower  layer  is  nearly  vertical.  One  deals  with  the 
departures  in  pressure,  temperature,  and  density  from  the  reference  state  given 
in  Sect.  2.2. 

In  the  general  mean  turbulent  departure  energy  equation  given,  the  presently 
relevant  term  is  the  turbulent  heat  flux  vector 

Q,  = p c . (2.3.1) 

3 o po  3 

In  the  general  departure  mechanical  energy  relation  for  the  mean  fluctuation 
turbulence,  the  presently  relevant  term  is  the  thermal  buoyancy  power  expended 
in  unit  volume 

p^gB^f^  . (2.3.2) 

The  lower  layer  is  gravitationally  unstable,  so  that,  in  the  atmosphere  and 
ocean,  the  Rayleigh  number  would  be  sufficiently  high  to  give  immediate  transi- 
tion to  turbulence.  The  buoyant  parcels  of  fluid  shown  in  Fig.  2 are  turbulent 
eddies;  whether  the  eddies  move  up  or  down,  the  associated  turbulent  heat 
correlation 


9 


T*u^  > 0 


for  unstably  stratified  fluid,  whereas 


T'u'  = 0 


(2.3.4) 


to  excellent  approximation  for  the  stably  stratified  fluid  flow.  In  Fig.  2 it 
might  be  supposed  that  a loose  turbulent  analogue  to  the  Fourier  heat  law 
prevails  so  that  ~ k^3T3x^ , where  is  a Boussinesq  eddy  conductivity. 

But  this  loose  analogue  ignores  the  connection  between  buoyancy  force  and  heat 
transfer.  Moreover,  the  turbulent  heat  flux  correlation  T'uj  does  not  depend 
on  the  local  T and  its  derivatives  (tracking)  but  has  an  identity  of  its  own 
in  the  fluid. 

Let  us  examine  the  two  implications  of  the  preceding  paragraph  in  Fig.  2. 

In  stratum  P,  the  pause  layer,  we  see  that  it  is  heated  from  below  but  no  heat 

leaves  it  from  above.  The  turbulent  heat  correlation  has  switched  off.  We  then 

ask,  what  is  the  fate  of  the  heated  pause  layer?  In  the  extensive  lower  unstable 

layer,  the  turbulent  heat  flux  T*u'  is  clearly  driven  by  f(x.)  in  a bulk 

o o 

fashion  but  there  is  no  one-to-one  correspondence,  in  general,  between  the  two. 
The  stratum  L shown  in  Fig.  2 is  heated  from  below  by  the  turbulent  heat 
correlation  and  cooled  from  above  thereby,  but  not  to  the  same  extent;  we  cannot 
expect  that  3(T'u^)/3Xj  = 0 in  general.  In  both  the  strata  P and  L we  see 
a turbulent  heat  correlation  imbalance  that  must  be  resolved. 

To  this  end,  we  turn  to  the  general  mean  turbulent  energy  equation.  If 
we  neglect  therein  the  viscous  dissipation,  the  viscous  turbulent  dissipation, 
the  molecular  conductivity,  and  radiative  transfer,  we  arrive  at  the  buoyancy 
equation  (that  governs  the  scalar  part  of  the  departure  buoyancy  force  per  unit 
volume  in  the  momentum  equation),  as  follows: 

h ^ I?  * sir  ^ <3  ' 

(2.3.5) 


For  the  initial-value  problem  under  discussion  in  Fig.  2,  we  have  first  to  look 
at  the  term  corresponding  to  turbulent  heat  flux  in  the  j = 3 direction,  viz. 


(2.3.6) 


Evidently,  this  term  is  positively  infinite  in  the  pause  layer  (if  the  pause 
layer  is  thin),  nonzero  in  the  lower  layer,  and  zero  in  the  upper  layer.  On 
the  other  hand,  the  convective  term  in  (2.3.5)  is  satisfied  everywhere  by  the 
initial  conditions  given  in  Fig.  2.  The  nonvanishing  of  term  (2.3.6)  in  the 
pause  layer  and  the  lower  layer  is  crucial;  we  see  in  (2.3.5)  that  it  leads 
initially  to  local  heating  (or  cooling)  because  we  must  get 


3T/8t  i 0. 


In  the  momentum  equation,  if  the  departure  buoyancy  terms  'v  T becomes  fully 

time  dependent,  then  hydrostatic  equilibrium  in  the  x direction  becomes 

o 

impossible  so  that  mean  motion  ensues  in  the  lower  layer. 

To  see  what  sort  of  mean  motion  results  in  the  steady-state  operation, 
suppose  that  the  flow  is  approximately  steady  in  a moving  frame  attached  to 
Then  (2.3.5)  gives 

(p  6 gfu.  + p 6 gf^)  = 0.  (2.3.7) 

3x.  o o 1 o o 1 
3 

This  conservation  equation  means  that  the  nonvanishing  of  (2.3.6)  must  be 
counterbalanced  by  considering  other  components  of  the  correlation  to  be 
nonzero  and  especially  by  mean  convective  motion  in  the  interior  of  the  lower 
layer. 

__  0 

Suppose  u,  = T'u'  = T — = 0 at  first.  Then  j = 3 in  (2.3.7)  and 
11  3x^ 

(2.3.7)  is  integrable,  to  give 

T(x2)u3(x^)  + T'u^  (Xg)  = T'u^  (Xg  = 0), 
if  u (x  ) = 0.  The  result  signifies  that  mean  convective  motion  extends  from 

O sj 

the  pause  layer  vertically  downward  to  infinity.  Clearly,  this  result  is  crude 
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and  we  should  consider  next  the  case  where  u^,  and  3/3x^  are  nonzero. 

If  we  examine  (2.3.5),  using  the  substantial  derivative  on  the  left  and  the 
positive  value  of  term  (2.3.6)  on  the  right,  then  at  the  top  of  the  lower  layer 
near  the  pause  layer  we  see  that  T for  a moving  fluid  element  must  increase. 

The  two  considerations  of  the  last  paragraph  lead  us  from  vertical  con- 
vection to  the  type  of  cellular  convection  sketched  in  Fig.  2,  where  the  cells 
corotate,  not  counterrotate.  This  cellular  motion  in  the  extensive  lower  layer 
is  the  direct  consequence  of  the  nonvanishing  of  term  (2.3.6)  in  the  pause  layer 
and  lower  layer.  Noting  that  we  deal  always  with  departures  from  the  reference 
state,  we  observe  that  the  cellular  motion  implies  nonzero  values  of  T'u|  and 
T'u'  as  well  as  T'u'  > 0. 

/.  <>3 

We  may  now  generalize  Eqs.  (2.3.3)  and  (2.3.4)  to  read  that,  for  the  heat 
switch, 

T'ul  > 0,  T'u'  0,  (a  = 1,2)  (2.3.8) 

3 a 

for  unstably  stratified  flow,  whereas 

TMlJ  = 0,  (j  = 1,2,3)  (2.3.9) 

to  excellent  approximation  for  the  stably  stratified  flow. 

2.4  Heat  Engine 

At  the  structured  interface  (called  the  pause  layer)  in  Fig.  3 two  events 
happen.  First,  the  thermal  cells  and  turbulent  eddies  buffet  the  stable  upper 
layer  by  "convective  penetration"  (the  phenomenon  was  first  noted  in  a tank 
of  water  by  Townsend  1964),  delivering  power  to  the  upper  layer  by  the  process. 

Moreover,  it  is  to  be  borne  in  mind  that  the  turbulent  heat  flux  must  switch 
off  rather  abruptly  upon  reaching  strong,  stable  stratification.  The  second  event 
at  the  interface  is  therefore  the  formation  of  a thermal  layer  that  can  lower  the 
local  Richardson  number  in  a thin  layer  from  a stable  value  to  a value  less 

than  the  critical  value  0.25  for  which  the  layer  becomes  dynamically  unstable. 

j 


The  two  events  produce  respectively  internal  waves  and  Kelvin-Helmholtz  waves. 
The  breaking  of  these  waves  ultimately  absorbs  the  thermal  buoyancy  power 
delivered  at  the  interface  from  below.  These  two  phenomena  are  depicted  in 
Fig.  3.  Work  is  therefore  being  done  on  the  upper  stable  layer  by  the  thermal 
phenomena  in  the  lower  layer  so  we  have  what  might  be  called  a "heat  engine." 

3.  The  Specific  Problem 

Our  prior  experimental  probing  in  a low-speed  wind  tunnel  (cf.  Appendix) 
with  our  concomitant  theoretical  considerations  suggest  the  common  existence  of 
a system  of  large  and  small  parameters.  These  plus  other  considerations  permit 
the  formulation  of  our  relatively  simple,  rational  theory  of  thermally  turbulent 
incident  stratified  flow  that  is  turbulent  due  to  gravitational  instability. 
Specifically,  outside  any  interfacial  pause  layer,  turbulent  stresslike  terms, 
turbulent  kinetic  energy  terms,  and  viscous  dissipation  terms  can  be  removed 
mainly  because  of  certain  large  correlation  ratio  parameters  Jp,  I,  and  H 
to  be  defined  later;  moreover,  it  is  common  that  the  relevant  values  of  the 
reciprocal  internal  Froude  number  are  small  and  it  is  convenient  to  our  purposes 
to  develop  the  square  of  the  reciprocal  internal  Froude  number  in  ascending 
series  in  this  section  as  an  added  simplification  and  elucidation.  We  present 
in  Sect.  3 the  theory  and  its  application  to  the  problem  that  will  now  be 
described. 

3.1  Physical  Statement  of  Specific  Problem 

We  return  now  to  Fig.  1 and  the  discussion  in  Sect.  2.1.  We  undertake 
the  considerations  of  turbulent  flow  due  to  large  P (but  not  large  R ) 

3 6 

due  to  the  gravitationally  unstable  region  depicted  therein;  the  equations  also 
hold,  often  trivially,  in  the  nonturbulent  stable  region  where  the  fluctuations 
are  identically  zero  so  that  the  entire  fluid  domain  is  treated  except  the 
interface  in  Fig,  1 when  the  thin  pause  layer  structure  of  the  interface  must 
be  considered  (see  Fig.  2). 
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For  X = - d and  for  x = + D,  we  provide  boundary  conditions  of  zero 
heat  transfer  on  the  solid  boundaries  over  the  interval  0 < x^^  < ».  At  x^^  = 0, 
the  depicted  temperature  profiles  T(Xg)  and  velocity  profiles  'the 

upstream  boundary  conditions  introduced,  superposed  on  which,  and  not  shown,  are 
suitable  turbulence  fields.  For  x^  = 0,  the  boundary  condition  is  that  turbu- 
lence vanishes  in  the  stable  layer.  In  the  unstable  layer,  the  boundary  condition 
at  x^  = 0 is  taken  to  be  that  of  turbulent  thermal  convection  produced  by  an 
extensive  uniformly  heated  surface  for  small  Prandtl  number  P^,  large  R^, 
large  flux  Richardson  number  J^,  large  I,  and  large  H.  Of  course,  no 
requirement  is  set  on  internal  Froude  number.  The  requisite  turbulent  correlations 
can  be  taken  from  experiment  and/or  theory.  These  data,  in  truth,  fix  the  problem. 
The  boundary  data  are  then  set  in  semi-infinite  parabolic  form.  The  problem  has 
its  antecedents  in  the  highly  nonlinear  extension  of  the  B^nard  problem  but  is 
aperiodic  in  x^,  driven  only  by  the  upstream  boundary  conditions,  and  approaches 
some  equilibrium  solution  as  x^^ 

3.2  Formulation  and  Procedures  of  Specific  Problem 

3.2.1  Ansatz 

We  begin  with  the  general  unsteady  laminar  departure  Eqs.  (2.2.5),  (2.2.6), 
and  (2.2.7)  with  p =0,  k =0,  F.  = 0,  and  (p  /2)e..e..  = 0.  The  energy 
equation  E reduces  to  the  "buoyancy  equation"  B that  amounts  to  the  substantial 
derivative  of  the  buoyancy  force  magnitude  in  M. 

The  departure  variables  are  u^,  f,  and  f and  the  Boussinesq  thermal 
state  equation  remains  as 

p = p (1  - B T)  (3.2.1) 

o o 

(For  a perfect  gas,  3 = 1/T  .) 

o o 


j 

J 

i 

li 


i i 


We  decompose  the  departure  variables  into  the  two  parts: 


k 


TCx,) 

+ 

T(x^,  x^,  Xg,t) 

ntxj) 

+ 

p(x^,  X^,  Xg,t) 

+ 

Uj(Xi,X2,X3.t) 

Sheared 

"Excursion  Flow" 

hydrosta- 

due to  cells. 

tic  given 

turbulence,  and 

boundary 

conditions 

waves 

This  decomposition  serves  to  introduce  the  boundary  conditions  into  the  differen- 
tial equations  as  well  as  the  Brunt-VSisaiS  parameter 


0 3T  - 


We  form  dimensionless  equations  from  the  definitions 


U t 


u*  = _i.  p*  = _E — T*  = 2_  U*  = — , t*  = X*  = 

1 U ’ P „2’  ^ ^ U ’ d ’ 1 d 

c p U T c 

o c c 


(3.2*3) 


where  U is  a constant  characteristic  velocity  difference  in  the  x direction 


d is  a constant  characteristic  associated  distance  in  the  x direction 


and  T is  a constant  characteristic  value  of  T. 
c 


We  obtain  the  following  excursion-flow  equations  in  the  dimensionless  excur- 


sion-flow variables  u’|,  p*,  and  T*: 


C*: 


^ = 0, 
3 ^ 


(3.2.4) 


M*: 


H*  (0«  . U‘(x3)  ^ <Gp  . U..(x5)  . ^ (S«*) 

1 3 3 J 


= ; (i,j  = 1,2,3), 


(3.2.5) 


and 


B*: 


X*)/£lu*. 


J(x* 


(3.2.6) 
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For  similitude  parameters,  we  find  above 


-2  -2  - 
F _ 5 — = e = (internal  Froude  no.)  > 0, 

^ ^ U 

c 


(3.2.7) 


N^(x„) 

J(x.)  E ^-4- 

^ U^/d^ 

c 


a Richardson  no-.  - 0. 


(3.2.8) 


Now  is  small  in  the  Boussinesq  approximation,  tending  to  make  e 

small.  Moreover  |j/e|  1 for  B*  to  balance  for  small  e.  Hereafter  we 

shall  exploit  by  series  expansion  the  small  parameters 


c = 10 


and  I d/e  I 1» 


(3.2.9) 


noting  that  the  smallness  of  these  parameters  affects  the  magnitude  of  the  mean 


flow  variables  but  not  the  relevant  turbulent  correlations  = 0(e  ) introduced 


by  the  upstream  boundary  conditions  at  full  strength  for  all  e.  The  relative 
importance  of  the  turbulent  correlations  is  determined  by  additional  considerations 
for  certain  large  correlation  ratio  parameters  Jp,  1,  and  H. 

Henceforth,  let  the  notation  simplify  to 


(3.2.10) 


3.2.2  Mean  Turbulent  Equations 


We  decompose  the  excursion  flow  into  cells  or  waves  ( ) plus  temporal 


I 
I 


turbulent  fluctuations  ( )'  according  to 


u.  = u.  + u!,  P = P + p',  T = T + T', 

and  time  mean  the  preceding  C,  M,  and  B to  obtain 


(3.2.11) 


k,  k'  * 


u. , . = 0, 
: D 


B:  + U{x3)G.,^  + ^ ' ^’i’ 


(i.j  = 1,2,3),  (3.2.13) 


f,^  + UCXg)?,^  + (fGj,j  = - (T'up,j  - |j(x3)/ep3. 


(3.2.14) 


The  variables  are  u^,  p,  T,  and  (T'up  and  the  equations  do  not  close 
because  the  three  correlations  (T’u^)  are  as  yet  unspecified.  In  the  mean 
equations  of  motion,  a necessary  condition  for  neglecting  the  Reynolds'  stress 
term  in  M is  that  the  flux  Richardson  number  Jp  is  large. 

3.2.3  Series  Procedure  in  e for  Solving  Mean  Turbulent  Equations 
(a)  Zeroth-order  solution: 

If  € = 0,  then  C and  M uncouple  from  B.  Considering  the  boundary 


conditions,  the  solution  is 


u!°^  = 0 and  p^°^  = 0. 
1 


(3.2.15) 


As  long  as  turbulence  is  present,  B is  nontrivial  and  B becomes  the  following 
brief  equation  in  and  (T'up^°^: 


+ U(x  )fi°^  = - (T'uI)S?^  = 0(1) 

t j j.  3 3 


(3.2.16) 


on  measure  e.  This  relation  shows  no  mean  mechanical  activity  but  a mere  heating 
. 2 

up;  notice  that  J N is  absent  and  the  interface  can  merely  rise  and  erode 
the  upper  stably  stratified  layer. 

(b)  First-order  solution: 

If  0(e)  <<  1,  then  C and  M are  weakly  coupled  to  B and  B is 
modified  by  the  factor  pL  + 0(e)*].  It  will  be  observed  that  J "v  now  appears. 
Now  wn  arrive  at  the  peculiar  series 


■0-r 


I I 


I 


i 
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^ = 


1 


+ ... 


P = 
T = 


0 

;(o) 


f<»-c 


+ . . . 


(3.2.17) 


ir  • • • 


(T'up  = (T'u!)^°^  + (T'up^^^*e  + 


Then  the  first-order  equations  are 


= 0, 

3 3 


B‘^>:  G‘i>  t . G<^>«^.U,3(X3>  . 


(3.2.18) 

(3.2.19) 


and 


- B<°>; 


= - (T'upS^^  - [j(Xg)/E3u^^\ 


(3.2.20) 


The  first  two  equations  of  this  trio  uncouple  from  the  third;  they  are  linear 


and  we  solve  for  and  Subsequently,  we  are  able  to  solve  the  thira 


?(1) 


equation;  it  is  linear  and  we  ultimately  solve  for  and  (T'u^)^^^. 


3.2.4  Mean  Turbulent  Fluctuation-Moment  Equations 


Begin  with  C,  M,  and  B,  Eqs.  (3.2.4)  to  (3.2.6),  write  u.  = u,  + u'. 

= = = 111’ 


T = T + T',  etc.,  and  next  multiply  M and  B each  by  Uj^  and  then  T’ 

to  form  the  following  mean  turbulent  fluctuation-moment  equations  indicated 
in  symbolic  form  by 


M*Uj^,  M-T',  B*T’,  and  B«u^.  (3.2.21) 

In  closing  the  mean  turbulent  equations  for  solution,  we  work  not  with  the 
traditional  but  with  the  three  equations 

[y*T’  + puj]  (3.2.22) 

and  the  one  equation 


[B*T'] 


(3.2.23) 


1 


A 
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governing  as  variables  the  two  appearing  correlations  (T’u|)  and  (T’T’);  we 


-/ 


neglect  (uTuT)  compared  to  (T’up  as  they  occur  in  (3.2.22)  and  justify  this 
strongest  assumption  made  by  use  of  the  argument  mentioned  in  the  last  paragraph 
of  Sect.  1 and  attributed  to  Lumley  and  Khajeh-Nouri  (1973). 

For  the  "extra"  correlations  in  these  fluctuation -moment  equations,  we 
arbitrarily  adopt  from  the  literature  the  turbulence  modeling  of  Donaldson  (1973) 


for  (p'T'),  (p'T}^),  (T'u|Uj),  and  (T'T'up;  this  modeling  introduces  the 

length-scale  parameters  ^^(oi  = 1*2,3)  and  the  velocity  parameter  V.  We 
emphasize  that  our  contribution  in  this  paper  is  not  dependent  on  the  individual 
choice  made  or  preferred  for  the  turbulence  modeling.  Using  the  nondimensionalization 
method  of  Sect.  3.2.1,  the  above  cited  turbulence  modeling  gives 


(p'T')  = - A V(T'uI),. 

J J 


(p'T;^)  = - (V/Aj^)(T’up 


(3.2.24) 


(T'ujuj)  = - VA^J^d'up,^  + (T’ut),J 


(T'T'up  = -VA^(T'T'),^ 

The  four  fluctuation-moment  equations  discussed  above  are  now  written 


explicitly.  From  [M*T'  + B*u|]  we  obtain  the  three  explicit  equations 


(T'up,^  + u^(T'up,^  + U(x3)(T'u|),^  + 6^^U,3(x3)(T'u^)  + u^,j(T'uJ) 


- {A-VC(T'u'.),,  + (T'u'.),.]},.  = e6,.(T'T')  t [A,V(T 'u  1 ) , . ] , . - j-  (T’ul); 


1 '3 


3 1 3 


3i  ‘•  3 3 'j-'i 

(i,j  = 1,2,3), 


(3.2.25) 


and  from  [B*T']  we  get  one  explicit  equation 


f (T'T'),^  + J u^(T'T'),.  + J U(x3)(T'T'),^  + T,^(T'uj)  - ^ [A2V(T 'T ' ) , ^ ] , ^ 

(3.2.26) 


- [J(x3)/e](T’u^), 
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The  four  fluctuation -moment  equations  above  plus 
Q,  M,  B close  the  thermal  turbulence  system  under 
study.  The  unknowns  are  u^,  p,  T,  (T'u|),  and 
(T'T').  (3.2.27) 


3.2.5  Series  Procedure  in  e for  Solving  Mean  Turbulent  Fluctuation- 
Moment  Equations 
(a)  Zeroth-order  solution: 

If  e = 0,  we  have  already  found  that 


u^°^  = 0 and  p^°^  = 0. 
1 ^ 


(3.2.15  bis) 


Then  the  fluctuation-moment  equations  above  give  for  [M*T'  + B«u!]^°^  the 


result 


(T'uj)S°^  + U(x3)(T'upi°^  + 6^j,U,3(x3)(T'up^°^  - {A2V[(T’up,j  + (T'up,^])S°' 


= + [A.V(T'un,,]S°^  - [ (T^)]^°\  (3.2.28) 

o J J i A 


and  for  [B*T']^°^  the  result 


i(T'T’)S°^  + i U(x  )(T'T')S°^  + fi?\T’up^°^  - i [A-V(T'T'),.]S?^ 


= [J(x3)/e3(T'u')^°\ 


(3.2.29) 


The  first  of  this  pair  is  uncoupled  from  the  second.  The  first  is  to  be  solved 
simultaneously  with  Eq.  (3.2.16)  for  B^°^  for  the  variables  and 

(T'up^°\  After  this  solution,  the  second  of  this  pair  can  be  solved  for 
(t7tT)(o)  if  desired. 

(b)  First-order  solution: 


The  complete  series  is 


e 


>e  -i 

T = + ... 

(f^)  = + ... 

(YTfr)  = + (T^)^^^*e  t ... 


U.2.30) 


The  moment 


equations  give  for  [M*T'  + - [M*T*  + the  three 


= 1 


= 1 


following  linear  equations  in  the  variables  (T'u|)^^^  and 


(T'u')iJ;^  + G$^^t’u:)S°^  + u(x-)(T'u:)S?-^  + 6,  .u,,(x,)(t»u')^^^  + uj^^d’up^®^ 

iX  3 ^3  <s  o X j j 


{A-V[(T’u.'),.  + (T'u!),.]}Sj-^  = 6,.(T’T')^°^  + [A.Vd'uH , - [|-  (T'ul)]^^^ 

* 13  3^3  33^ 


1 

(3.2.31) 


In  the  above,  from  [B*T']"^''  under  (a),  (T'T')^°^  is  given  by  the  prior  result 


(o) 


1 (T'T')i°^  + I U(x3)(T'T')S°^  - I [A^VCT'TO.jlSj 


* - - [J(x3)/ej(T'u’)^°\ 


(3.2.29  bis) 


4.  Evaluation  of  Power  Delivered  to  Upper  Layer  by  Heat  Engine 


Before  considering  from  Sect, 3. 2,  the  solution  equations  we  can  give  some 
explicit  results  for  the  heat  switch  of  Sect.  2.3  and  the  heat  engine  of  Sect.  2.4 
as  they  pertain  to  the  second  paragraph  of  Sect.  2.1  and  to  Fig.  3. 

Let  us  write  out  the  first  system  of  fluctuation -moment  equations  given  in 
Eqs.  (3.2.21).  These  were  indicated  in  Sect.  3 as  part  of  the  dimensionless 
excursion-flow  discussion  but  not  displayed  or  used.  For  M*u^  then  we  have  the 
relevant  Reynolds'  stress  equations 


“W-t  ♦ * -“j  ♦ "j“rj“k 


P’i^k  ’ 


(4.1) 
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After  solution,  we  may  return  to  this  equation,  contract  k to  i,  and 
obtain  the  mechanical-energy  relation  for  the  mean  fluctuations  according  to 
our  dimensionless  excursion-flow  equations.  Suppose  we  neglect 


(4.2) 


on  the  grounds  that  they  are  presently  and  commonly  small  as  argued.  We  then 


obtain  without  any  modeling  of  the  correlation  P»j“j  result 

= e(TX). 

I D 3 

Now  integrate  over  the  finite  volume  x and  use  the  identity 


(4.3) 


to  obtain  the  result 


b, .dx  = o bn.dA 
3 D 


(p'ut)n.dA  = - e (T'u’)dx 
3 3 ^ 

A •'X 


(4.4) 


for  a finite  fluid  aggregate. 

Therefore,  the  mechanical  power  delivered  to  the  upper  stable  layer  at 
the  interface  is  equal  to  the  thermal  power  of  the  thermal  turbulence  in  the 
lower  layer  and  vanishes  if  e = 0 (a  result  self-consistant  with  our  other 
results) . 

This  result  may  be  compared  to  the  dimensional  general  mean  fluctuation 
departure  energy  relation  (2.2.11)  above.  By  comparison,  we  can  infer  that 
the  assumptions  of  our  excursion-flow  theory  are 


P^6^g(T*up 

(l/2)li.ei. 


-6^g(T’up 


D_  'X 
Dt  2 


» 1, 


» 1, 


(4.5) 


(4.6) 
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H = 


(p'Uj ) 


(p  /2)(u!u!u!) 
o 1 1 j 


» 1, 


(4.7) 


plus  the  assumptions  of  negligible  viscous  dissipation 


-p  e g(T'u!) 
o o 3 

(u  /2)(e.’  .e.' . ) 
o i:  13 


(4.8) 


and 


•(p'uj ) 


y (e I .u ! ) 
o 13  1 


» 1. 


(4.9) 


5.  The  Solution  Equations  with  Arbitrary  Shear  U = U(x„) 

o 

The  requisite  equations  were  developed  in  Sect.  3 but  are  not  yet  in  a 
convenient  form  for  manipulation.  It  is  noted,  however,  that  all  equations 
are  linear  — in  some  cases  with  variable  coefficients  — even  with  arbitrary 
shear . 

5.1  The  Zeroth-Order  Equations  in  e 

We  gather  up  the  results  of  Sects.  3.2.3(a)  and  3.2.5(a).  We  have  noted 

2 

that,  to  zeroth  order,  there  exists  no  J N , Let 

$ - = (f^IJ)^°\  0 = (5.1.1) 

and 

Aj^,  A^,  AgjV  = assumed  consts.  (5.1.2) 


The  last  four  are  the  empirical  functions  of  Donaldson  (1973). 
Then  from 


®’t  ^^’‘3^®’! 

or,  in  vector  notation, 


(5.1.3) 


1 


If  . UU>  If  = - V? 


(5.1.4) 
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From  [M*T'  + 


= 1 


Q..^  t U(x3)Q,,^  , «U«.3<>‘3)'33  = * <‘2  * 

(5.1.5) 

or,  in  vector  notation, 


1^  + U(z)  U + 1 U(z)k*Q  = A^VV^Q  + (A^  + Ag)W(7  • $)  - (V/A^)Q 


(5.1.6) 

Given  U(z),  Eq.  (5.1.6)  can  be  solved  as  it  stands  for  We  then  solve 

Eq.  (5.1.4)  for  the  other  unknovm,  viz  0. 

5.2  The  First-Order  Equations  in  e 


We  gather  up  the  results  of  Sects.  3.2.3(b)  and  3.2.5(b).  We  repeat  that, 


to  first  order,  there  appears  a J 'u  N'^.  Let 


q * q.  = (T'u!)'^\  T = (T'T')""% 


0 — 


:(1) 


, V - V.  = u. 

1 X 


. = 5<“. 


(5.2.1) 


and 


A^jA^jA^jV  = assumed  consts. 


(5.2.2) 


n(l) 


Then  the  problem  is  given  by  solving  simultaneously  the  equations  C , 
snd  [M*T'  + B’u!]^^^  - [M*T'  + B*u!]^°^  since  the  equations  B^°^] 


= 1 


i.T 


= 1 


and  [B*T']  are  to  be  solved  aside.  We  therefore  write 


v.,.=0  or  7*v=0, 

: ] 


Vi>t  t U(X3)V.,^  t V36^.U,3(x3)  = 063.  - ir,. 


(5.2.3) 

(5.2.4) 


or 


and 


(5.2.5) 


[y*T'  + B-uj]^^^  - [y*T'  + B*up^°^ 

'^i’t  " ^j'^i’j  ^ ^^^3^‘li’l  ^ ^ - V'ii’jj 

= 63.T  + (A^  + A3)Vq  , - (V/A^)q.  (5.2.6) 


I 


+ (v*V)Q  + U(z)  ^ U(z)(k*q)  + (Q  • V)v  - A^VV^q 

= kT  + (A^  + Ag)VV(V-q)  - (V/A^)q.  (5.2.7) 

The  equation  for  t can  be  solved  alone  once  the  right  member  has  been  determined 
from  the  zeroth-order  solution.  Then,  from  [B*T']^°\ 


-iV'-Jj 


- Q.0,,  - [J(x-)/e]Q_  (5.2.8) 

J J o o 


7 It  T |J  ~ 7 ^ ■ t'J(2)/Ok*Q. 


(5.2.9) 


The  variable  9 does  not  appear  in  the  above  equations  but  the  buoyancy  equation 


depends  on  v and  q obtained  above.  To  be  solved  last  then  is 

1^^^  - 6,^  t U(Xg)9,^  +0,^Vj  + [J(x3)/e]Vg  = - .j 


(5.2.10) 


+ U(z)  1^  + (70)‘V  + [J(z)/e]ic*v  = - V*q. 


Discussion 


(5.2.11) 


What  is  now  required  is  to  solve  the  equations  of  Sect.  5 subject  to  the 
boundary  conditions  given  plus  the  boundary  conditions  at  the  upper  lid  or 
interface  near  x^  = 0.  The  difficulty  is  that  the  upper  layer  and  lower  layer 
are  closely  coupled  together;  although  the  lower  layer  is  the  driver  as  per 
Sect . 4 , there  will  be  feedback  from  the  upper  layer . Thus  we  must  deal  above 
with  internal,  surface,  and  Kelvin-Helmholtz  waves  simultaneously.  Then  a wide 
variety  of  forms  of  cellular  structure  below  can  be  obtained  by  suitable  adjust- 
ment of  the  shear  rate  U(x  ).  A survey  of  the  limited  literature  on  this 
coupled  problem  is  given  by  Turner  (1973);  see  Ch.  7 and  especially  paragraph  7.3, 
p.  226.  There  is  an  extensive  literature  on  linear  waves  in  stably  stratified 
flows  and  wo  have  rendered  our  unstably  stratified  lower  turbulent  flow  linear 
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and  therefore  compatible.  (The  problem  is  still  well  formulated  if  the 
expansion  in  e is  not  made.)  The  detailed  solution  is  left  for  further 
research. 

This  research  was  supported  by  the  Fluid  Dynamics  Program  of  the 
U.  S.  Office  of  Naval  Research  under  contracts  N0001‘^'-67-A-0191-0021  and 
N00014-76-C-0154.  Thanks  are  due  to  S.  K.  F.  Karlsson  and  R.  S'.  Scotti  in 
connection  with  the  following  Appendix. 

Appendix:  Prior  Experiments 

We  have  mentioned  that  the  work  proposed  is  backed  up  by  prior  experi- 
mental probing  in  a low-speed  wind  tunnel  at  Brown.  The  experiments  were 
exploratory  and  thus  produced  qualitative  smoke  photographs  resembling  Fig.  3 
as  well  as  semi-quantitative  tunnel  data  such  as  shown  in  Table  1 to  be 
discussed  below.  The  research  was  reported  in  Clarke  et  al.  (1973)  and 

Clarke  (1975).  The  semi-quantitative  data  (only)  corresponded  physically  to 

-2  • 

the  case  e H F^  =0  in  the  sense  that  the  general  resemblance  was  remarkable; 

as  is  indicated  in  Eq.  (4*^),  the  heat  engine  did  not  operate  detectably;  the 

phenomena  observed  was  thermal  in  nature  with  little  mechanical  activity  of 

2 

wave  or  cell  type  (in  the  excursion  theory,  N is  absent  if  e = 0).  On 
the  other  hand,  the  qualitative  data  displayed  a great  deal  of  wave  activity 
of  all  types  and  suggested  strongly  the  operation  of  the  heat  engine  as  in 
Fig.  3. 

Both  the  qualitative  and  semi-quantitative  data  indicated  strongly  that 
the  flow  was  essentially  free  of  stress-like  turbulent  correlations  and 
dissipation  outside  of  the  thin  pause  layer  (and  wall  boundary  layers)  over 
a test  length  of  25  feet  within  a tunnel  cross-section  of  32  inches  in  width 
and  22  inches  in  height.  If  x^  = x in  the  flow  direction,  x^  = y in  the 
horizontal  direction  measured  from  one  vertical  wall,  x^  = z measured  from 
the  tunnel  floor,  and  since  6^^  = T^  for  a perfect  gas,  then  we  can  display 
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some  of  the  semi-quantitative  data  in  Table  1.  The  lower  unstable  layer 

corresponds  to  0 < z < 12,19  inches.  All  of  the  symbols  have  already  been 

“2  -2  -2 

defined.  It  is  seen  that  e = <10  , |J/Fj  | “v  1,  and  the  four  measured 

values  of  [j^l  are  comfortably  greater  than  unity  (the  values  not  shown  are 
larger).  All  of  the  hot-wire  data  taken  displayed  a decay  of  activity  at 
points  more  and  more  above  the  pause-layer  zone  of  activity. 

-2 

It  was  not  possible  to  raise  the  quantitative  values  of  e = Fj  to  the 
more  interesting  range  of  values,  clearly  perceived  theoretically,  from  the 
tunnel  results  given  above  in  the  Brown  University  wind  tunnel  without  sacrifice 
of  flow  uniformity.  On  the  other  hand,  the  experimental  and  theoretical  probing 
seemed  to  open  up  entire  new  theoretical  avenues  for  unstable  turbulent  flow  -- 
not  only  for  the  specific  problem  just  discussed,  but  for  similar  ones.  None  of 
this  probing  suggested  that  the  large  and  small  parameters  observed  were  limited 
in  occurrence  to  the  wind  tunnel  experiences  but  should  be  quite  common  in 
nature.  Accordingly,  we  undertook  to  develop  the  new  simplified  theory  of 
phenomenological  thermal  turbulence  suitable  whenever  the  parameters  behave  as 
required  and  subject  to  systematic  (i.e.  rational)  improvement  in  this  respect. 
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